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Abstract. Systematic uncertainties in the recent precise determination of the strangeness magnetic moment 
of the nucleon are identified and quantified. In summary, G S M = —0.046 ± 0.019 /ijv. 
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1 Introduction 

Recent low-mass lattice QCD simulation results combined 
with new chiral extrapolation techniques and the principle 
of charge symmetry have enabled a precise determination 
of the strangeness magnetic moment of the nucleon G S M 
1 . In this paper, the systematic errors of the approach are 
explored and quantified. In particular, we examine the sen- 
sitivity of G S M and the magnetic moments of the baryon 
octet to the regulator-mass scale of finite-range regular- 
ized chiral effective field theory, the lattice scale determi- 
nation, the finite-volume of the lattice, and the quenched 
approximation. 




u,d,s 



Fig. 1. Diagrams illustrating the two topologically different 
insertions of the current within the framework of lattice QCD. 
These skeleton diagrams for the connected (left) and discon- 
nected (right) current insertions are dressed by an arbitrary 
number of gluons and quark loops. 



2 Charge Symmetry 

The approach 1 centres around two equations for the 
strangeness magnetic moment of the nucleon, G S M , ob- 
tained from charge symmetry 
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Here the baryon labels represent the experimentally mea- 
sured baryon magnetic moments and R 8 d = G s M j l G d M is 
the ratio of s- and d-quark sea-quark loop contributions, 
depicted in the right-hand diagram of Fig. ^ l R S d nes m 
the range (0,1). The ratios u p /u E and u n /u~ are ratios of 
valence-quark contributions to baryon magnetic moments 
in full QCD as depicted in the left-hand diagram of Fig. Q 
The latter are determined by lattice QCD calculations ^ 



Equating Eqs. Q and @ provides a linear relationship 
between u p ' jvF and u n /u a which must be obeyed within 
QCD under the assumption of charge symmetry - itself 
typically satisfied at the level of 1% or better [5]. There are 
no other systematic uncertainties associated with this con- 
straint. Figure |21 displays this relationship by the dashed 
and solid line. Since this line does not pass through the 
point (1.0, 1.0), corresponding to the simple quark model 
assumption of universality, there must be an environment 
effect exceeding 12% in both ratios or approaching 20% or 
more in at least one of the ratios. To determine the sign of 
G S M , it is sufficient to determine where on this constraint 
curve, QCD resides. 
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Fig. 2. The straight line (dashed G S M (0) < 0, solid G S M (0) > 0) 
indicates the constraint on the ratios u p /u s and u n /u" implied 
by charge symmetry and the experimentally measured mag- 
netic moments. The assumption of environment independent 
quark moments is indicated by the square. The dependence of 
the extrapolated ratios from lattice QCD simulations on the 
parameter A — 0.7, 0.8 and 0.9 GeV, (governing the size of 
pion-cloud corrections) is illustrated by the cluster of points 
with A increasing from left to right. 



3 u p /u and u n /u~ Determination 

Our present precise analysis has been made possible by a 
significant breakthrough in the regularization of the chiral- 
loop contributions to hadron observables [7||^Hi?j . Through 
the process of regulating loop integrals via a finite-range 
regulator (FRR) 8,10], the chiral expansion is effectively 
re-summed to produce an expansion with vastly improved 
convergence properties. 

The chiral expansion for the w-quark contribution to 
the proton magnetic moment in quenched QCD (QQCD), 
has the form 



a + uP Xv' V ("V J A ) + XttB Ib (Wtt , A) 

+Xkb lB{m K , A) + «2 ml + of m%, 



(3) 

where the repeated index, B, sums over allowed baryon 
octet and decuplet intermediate states. The dependence 
of the unrenormalized coefficients, of, and the associated 
dipole- vertex regulated loop integrals, /(m W) yl), on the 
regulator parameter, A, is emphasized by the explicit ap- 
pearance of A. The loop integrals are denned as 
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where Abn is the relevant baryon mass splitting and the 
function u(k, A) is the dipole-vertex regulator. The coef- 
ficients, x, denote the known model-independent coeffi- 
cients of the non-analytic terms for tt and K mesons in 
QQCD [3111121- 



Figure illustrates a fit of FRR, quenched chiral per- 
turbation theory (%PT) to our lattice results (solid curve), 
where only the discrete momenta allowed in the finite vol- 
ume of the lattice are summed in evaluating the chiral loop 
integrals. The long-dashed curve that also runs through 
the lattice results corresponds to replacing the discrete 
momentum sum by the infinite- volume, continuous mo- 
mentum integral. For all but the lightest quark mass, finite 
volume effects are negligible. 

Having determined the analytic coefficients Oq 2 4 f° r a 
particular choice of A, one can correct the chiral properties 
of the pion-cloud contribution from QQCD to full QCD 9 

by changing the coefficients of the loop integrals, XV, 
XttB, Xkb of Eq. Q, to their full QCD counter parts ^2 
\F2\ . Valence quark contributions in full QCD are indicated 
by the long-dash-dot curve in Fig. (i.e. sea-quark loop 
charges are zero) and the full it-quark sector including 
the u-sea-quark loop contributions are indicated by the 
short-dash-dot curve for A = 0.8 GeV. Figures 0] [S] and 
|6|show similar results for the u quark in n, E + , and S° 
respectively. 

The importance of correcting for both finite-volume 
and quenching artifacts is illustrated in Figs. [3 and |H1 
where the one standard deviation agreement between the 
chirally corrected lattice QCD simulation results and the 
experimentally measured baryon magnetic moments is high- 
lighted. 



4 Systematic Errors 
4.1 Regulator Dependence 

It is important to investigate systematic errors associated 
with the regulator-mass dependence of FRR x?T • The 
extrapolated results of finite-volume quenched chiral ef- 
fective field theory should be insensitive to the choice of 
regulator parameter. When working to sufficient order in 
the chiral expansion, changes in the regulation of loop in- 
tegrals should be absorbed by changes in the unrenormal- 
ized coefficients, flQ 2 4, in a manner which preserves the 
invariant renormalized coefficients. The latter are reflected 
in Fig. which illustrates the insensitivity of quenched 
baryon magnetic moments in a finite-volume to the regu- 
lator parameter A. This systematic error is small relative 
to the statistical error. 

Since the finite-volume and quenching corrections are 
applied only to the loop integral contributions, the final 
results are A dependent. In this case, the regulator of 
the loop integral has become a model for the axial-vector 
form factor of the nucleon, describing the coupling of pi- 
ons to a core described by the analytic terms of the FRR 
expansion. This approach describes the relation between 
quenched and full QCD N and A mass simulation results 
as a function of m, very accurately [E]. FiguresITT)! andlTTl 
display results for positive and negative baryon magnetic 
moments in full QCD respectively. One-standard devia- 
tion agreement is achieved for 0.7 < A < 0.9. 

Because the strangeness magnetic moment of Eqs. {Q 
and J2Jl depends only on ratios of magnetic moments, most 
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Fig. 3. The contribution of a single u quark (with unit charge) 
to the magnetic moment of the proton. Lattice simulation re- 
sults (square symbols for > 0.05 GeV 2 ) are extrapolated 
to the physical point (vertical dashed line) in finite-volume 
QQCD as well as infinite volume QQCD. Estimates of the va- 
lence u quark contribution in full QCD and the full u-quark 
sector contribution in full QCD are also illustrated. Extrapo- 
lated values at the physical pion mass (vertical dashed line), 
are offset for clarity. 
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Fig. 5. The contribution of a single u quark (with unit charge) 
to the magnetic moment of S + . Curves and symbols are as for 
Fig El 
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Fig. 4. The contribution of the u quark (with unit charge) to 
the magnetic moment of the neutron. Curves and symbols are 
as described in Fig. El 



Fig. 6. The contribution of the u quark (with unit charge) to 
the magnetic moment of the S° hyperon. Curves and symbols 
are as described in Fig. El 



of this A dependence cancels in the final ratios, as illus- 
trated by the close clustering of points in Fig. 

4.2 Role of the Decuplet in %PT 

It is often argued that next-to-leading-order non-analytic 
(NLNA) contributions from the baryon decuplet are es- 
sential in describing the mass dependence of nucleon mag- 
netic moments. While the decuplet baryon contributions 
are not necessarily small, we find the non-analytic curva- 
ture induced by these contributions is sufficiently subtle 



that it may be accurately absorbed by the analytic terms 
of the chiral expansion. FigurelT^illustrates the insensitiv- 
ity of finite-volume quenched chiral effective field theory 
to NLNA decuplet-baryon contributions. 

Figure ^] confirms that the NLNA decuplet-baryon 
contributions are indeed large in some cases and as such 
are important in correcting the artifacts of the quenched 
approximation. However, other more highly-excited baryon 
resonances have small couplings to the ground-state baryon 
octet relative to that for the decuplet and provide negli- 
gible corrections. 
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Fig. 7. The one standard deviation agreement between the 
chirally corrected lattice QCD simulation results (square sym- 
bols) and the experimentally measured baryon magnetic mo- 
ments (circular symbols) having positive values. Finite- volume 
quenched results (crossed boxes) and infinite-volume quenched 
results (diamonds) are also illustrated to highlight the impor- 
tance of correcting for both finite-volume and quenching arti- 
facts. 



Fig. 9. The FRR x?T regulator mass dependence of finite- 
volume quenched chiral effective field theory. Experimental 
measurements are illustrated at left for each baryon for ref- 
erence. Results for A = 0.6, 0.7, 0.8, 0.9 and 1.0 GeV are illus- 
trated from left to right for each baryon. The small systematic 
dependence on A relative to the statistical error bars illustrated 
indicate the order of the chiral expansion is adequate for this 
analysis. 
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Fig. 8. The one standard deviation agreement between the 
chirally corrected lattice QCD simulation results and the ex- 
perimentally measured baryon magnetic moments having neg- 
ative values. Symbols are as in Fig. 



ticular, the Sommer parameter, ro, is ideal as it sets the 
scale by equating the force between two static quarks in 
QQCD and full QCD at a precise separation of r = 0.49 
fm. 

The string tension might also be used, but screening 
of the potential in full QCD makes this measure poorly 
defined. However, we consider it here as a measure of the 
systematic error encountered in setting the scale of the 
lattice QCD results. FigureH4l illustrates this scale depen- 
dence on baryon magnetic moments and Fig. ll5l illustrates 
the rather minor impact this systematic uncertainty has 
on the valence-quark moment ratios vital to determining 
the sign of G S M . 

Collecting the variation of the valence-quark moment 
ratios from variations in A (which maintain one-standard- 
deviation agreement with experiment), variations in set- 
ting the lattice scale and statistical errors determined by a 
third-order single-elimination jackknifc analysis, one finds 



4.3 Scale Dependence 



= 1.092 ±0.030 and 



= 1.254 ±0.124. 



(G) 



Setting the scale in quenched QCD simulations is some- 
what problematic. Different observables lead to different 
lattice spacings, a. If one is explicitly correcting the one- 
loop pion-cloud contributions to hadronic observables, as 
we are here, then clearly one must set the scale using an 
observable insensitive to chiral physics. This excludes ob- 
servables such as the rho-meson mass, nucleon mass, or 
the pion decay constant commonly used in the literature 
to hide the artifacts of the quenched approximation. On 
the other hand, the heavy-quark phenomenology of the 
static-quark potential provides an optimal case. In par- 



Using the experimental magnetic moments, one ob- 
serves that Eq. J2J) 



G\ 



e ns 

Kj 



-1.033 



0.599) 



(7) 



is least sensitive to variation in the valence-quark moment 
ratio, and hence provides the most precise determination 
for G S M . Figure HUI plots G S M as a function of l R s d with 
standard error limits provided by Eq. 1^. 
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Fig. 10. The dependence of positive octet-baryon magnetic 
moments on the parameter A — 0.6, 0.7, 0.8, 0.9 and 1.0 GeV, 
governing the size of pion-cloud corrections associated with 
the finite-volume of the lattice and artifacts of the quenched 
approximation. Experimental measurements, illustrated at left 
by the filled circle for each baryon, indicate that optimal cor- 
rections are obtained for 0.6 < A < 0.9 GeV. 
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Fig. 12. The insensitivity of finite-volume quenched chiral ef- 
fective field theory to NLNA decuplet-baryon contributions. 
Results including decuplet intermediate states in chiral effec- 
tive field theory (squares) are compared with results exclud- 
ing the decuplet (diamonds) . Experimental measurements (cir- 
cles), are illustrated at left for each baryon for reference. 
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Fig. 11. The dependence of negative octet-baryon magnetic 
moments on the parameter A = 0.6, 0.7, 0.8, 0.9 and 1.0 GeV, 
governing the size of pion-cloud corrections associated with 
the finite-volume of the lattice and artifacts of the quenched 
approximation. Experimental measurements, illustrated at left 
by the filled circle for each baryon, indicate that optimal cor- 
rections are obtained for 0.7 < A < 0.9 GeV. 



Fig. 13. The dependence of octet-baryon magnetic moments 
on the inclusion of NLNA decuplet-baryon contributions in the 
process of correcting the pion cloud of quenched chiral effective 
field theory. Preferred results including decuplet intermediate 
states in chiral effective field theory (squares) are compared 
with results excluding the decuplet (diamonds). Experimental 
measurements (circles), are illustrated at left for each baryon 
for reference. 



5 Estimating e R s d = G s M / e G d M 

The symmetry of the three-point correlation function |14j 
describing the sea-quark loop contributions to the nucleon, 
depicted in the right-hand side of Fig. ^ ensure that the 
chiral expansion for this quantity is identical for all three 
quark flavours, up to simple charge factors. For the d- or 
s-quark loop contributions, the only difference that can 
arise is whether one evaluates the chiral expansion at the 
pion or kaon mass. The leading non-analytic contribution 



to the chiral expansion involves two pseudoscalar meson 
propagators, and therefore one expects contributions to 
l R s d in the ratio m^/m 4 ^ ~ 0.1. 

To be more precise, one can use the same successful 
(single-parameter) model, previously used to correct the 
quenched simulation results to full QCD, as highlighted 
in Figs. [7| and 03 to provide an estimate for R s d - Evalu- 
ating the loop integrals with A = 0.8 ± 0.2 GeV yields 
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Table 1. Sources of uncertainty and their contribution to the strangeness magnetic moment of the nucleon, G 3 M , in units of 
nuclear magnetons, fj,N- Uncertainties are documented for G S M obtained from the valence-quark ratio u p /u s in Eq. Q, from 
the valence-quark ratio u n /u~ in Eq. (|2J and from a statistically weighted (SW) average of these two determinations. 



Uncertainty Source 


Parameter Range 


u p /u E , Eq. © 
G S M = -0.045 


U n /u 8 , Eq. 
G S M = -0.046 


SW Average 
G S M = -0.046 


Statistical Errors 




0.016 


0.009 


0.008 


Chiral corrections 


0.7 < A < 0.9 GeV 


0.001 


0.002 


0.002 


Scale Determination 


0.122 < a < 0.134 fm 


0.001 


0.002 


0.002 


i R s d Determination 


0.096 < 1 R s d < 0.181 


0.016 


0.017 


0.017 


Total Uncertainty 




0.023 


0.019 


0.019 
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Fig. 14. 


The dependence of octet-baryon magnetic moments 



on the matching criteria for determining the lattice spacing, a. 
While a r „ = 0.128 fm set by ro (square symbols) is the pre- 
ferred method for determining the scale, results for a CT = 0.134 
fm set by the string tension (triangles) and a third spacing of 
0.122 fm (diamonds) are also illustrated. Experimental mea- 
surements (filled circles) , are illustrated at left for each baryon. 

e R s d = 0.139 with 0.096 < e R s d < 0.181. This uncertainty 
dominates the final uncertainty in G S M . 

Table ^ summarizes the sources of uncertainty and 
their contributions to the final determination 

G S M = -0.046 ± 0.019 /j, N , (8) 

for the strange quark contribution to the magnetic mo- 
ment of the nucleon. 



Acknowledgments 

We thank the Australian Partnership for Advanced Com- 
puting (APAC) for generous grants of supercomputer time 
which have enabled this project. Support from the South 
Australian Partnership for Advanced Computing (SAPAC) 
and the National Facility for Lattice Gauge Theory is 
also gratefully acknowledged. DBL thanks Jefferson Lab 



1.3 



1.2 - 




0.8 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 2.5 

u n / u~ 

Fig. 15. The charge symmetry constraint line (dashed 
Gm(0) < 0, solid Gm(0) > 0) on the ratios u p /u s and u n /u s . 
The dependence of the ratios from chirally-corrected quenched 
lattice QCD on the scale parameter, a — 0.134, 0.128, and 
0.122 fm, is illustrated by the cluster of points with a decreas- 
ing from left to right. 



0.0 


i i i i i i i i 


-0.1 




-0.2 




^ -0.3 
— -0.4 




n 

^-0.5 




-0.6 




-0.7 




-0.8 


i i i i i i \i \ '.i i 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

a 

Fig. 16. The dependence of G S M on the strange to light sea- 
quark loop ratio l R a d . Standard error limits have there origin 
in the systematic error summary of Eq. ||3 



D. B. Leinweber et al: Systematic uncertainties in the precise determination. 



7 



for their kind hospitality where the majority of this re- 
search was performed. This work is supported by the Aus- 
tralian Research Council and by DOE contract DE-AC05- 
84ER40150 under which SURA operates Jefferson Lab. 



References 

1. D. B. Leinweber et al, arXiv:hep-lat/0406002 

2. G.A. Miller, B.M. Nefkens and I. Slaus, Phys. Rept. 194 
(1990) 1. 

3. J. M. Zanotti et al, Phys. Rev. D 65, 074507 (2002) 
|hep-lat/0U0216| . 

4. D. B. Leinweber et al, Eur. Phys. J. A 18, 247 (2003) 
|nucl-th/02lToi"4l . 

5. J. M. Zanotti, et a/., |hep-lat/0405015| 

6. J. M. Zanotti, et g/., |hep-lat /0405026 

7. R. D. Young et al , Prog. Part. Nucl. Phys. 50, 399 (2003) 
hep-lat/0212031 . 

8. D. B. Leinweber, A. W. Thomas and R. D. Young, Phys. 
Rev. Lett. 92, 242002 (2004) [hep-lat/0302020]. 

9. R. D. Young, D. B. Leinweber and A. W. Thomas, 
arXiv:hep-lat/0406001 

10. J. F. Donoghue, B. R. Holstein and B. Borasoy, Phys. Rev. 
D59, 036002 (1999). 

11. D. B. Leinwebe r, Phys. Rev. D 69, 014005 (2004) 
|hep-lat/0211017| . 

12. M. J. Savage, Nucl. Phys. A 700, 359 (2002) 



13. R. D. Young, D. B. Leinweber, A. W. Thomas 
and S. V. Wright, P hys. Rev. D 66 (2002) 094507 
|arXiv:hep -lat/0205017 . 

14. D.B. Leinwebe r, Phys. Rev. D53, 5115 (1996) 
|hep-ph/95123T9l 




